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O 1 Abstract 

] We construct self-similar functions and linear operators to deduce a self-similar variant of the Laplacian 
■ operator and of the D'Alembertian wave operator. The exigence of self-similarity as a symmetry property 
. . ', requires the introduction of non-local particle-particle interactions. We derive a self-similar linear wave oper- 
. ^ ■ ator describing the dynamics of a quasi-continuous linear chain of infinite length with a spatially self-similar 
^ distribution of nonlocal inter-particle springs. The self-similarity of the nonlocal harmonic particle-particle 
^ [ interactions results in a dispersion relation of the form of a Weierstrass-Mandelbrot function which exhibits 
- - ' self-similar and fractal features. We also derive a continuum approximation which relates the self-similar 
Laplacian to fractional integrals and yields in the low-frequency regime a power law frequency-dependence of 
the oscillator density. 
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2 Introduction 



In the seventies of the last century the development of the Fractal Geometry hj Mandelbrot [1] launched 
a scientific revolution. However, the mathematical roots of this discipline originate much earlier in 
the 19*'* century [2]. The superior electromagnetic properties of "fractal antenna^' have been known 
already for a while [3, 4]. More recently one found by means of numerical simulations that fractal 
gaskets such as the Sierpinski gasket reveal interesting vibrational properties [5]. Meanwhile physical 
problems in fractal and self-similar structures or media become more and more a subject of interest 
also in analytical mechanics and engineering science. This is true in statics and dynamics. However 
technological exploitations of effects based on self-similarity and "fractality" arc still very limited 
due to a lack of fundamental understanding of the role of the self-similar symmetry. An improved 
understanding could raise an enormous new field for basic research and applications in a wide range 
of disciplines including fluid mechanics and the mechanics of granular media and solids. Some initial 
steps have been performed (see papers [5, 6, 7, 8, 9, 10] and the references therein). However a 
generally accepted "fractal mechanics" has yet to be developed. Therefore, it is highly desirable to 
develop sufficiently simple models which are on the one hand accessible to a mathematical-analytical 
framework and on the other hand which capture the essential features imposed by self-similar scale 
invariant symmetry. The goal of this demonstration is to develop such a model. 

Several significant contributions of fractal and self-similar chains and lattices have been presented 
[13, 14, 15, 16]. In these papers problems on discrete lattices with fractal features are addressed. 
Closed form solutions for the dynamic Green's function and the vibrational spectrum of a linear chain 
with spatially exponential properties are developed in a recent paper [11]. A similar fractal type of 
linear chain as in the present paper has been considered very recently by Tarasov [7]. Unlike in the 
present paper the chain considered by Tarasov in [7] is discrete, i.e. there remains a characteristic 
length scale which is given by the next-neighbor distance of the particles. 

In contrast to all these works we analyze in the present paper vibrational properties in a quasi- 
continuous linear chain with (in the self-similar limiting case) infinitesimal lattice spacing with a 
non-local spatially self-similar distribution of power-law-scaled harmonic inter-particle interactions 
(springs). In this way we avoid the appearance of a characteristic length scale in our chain model. It 
seems there are analogue situations in turbulence [17] and other areas where the present interdisci- 
plinary approach could be useful. 

Our demonstration is organized as follows: § 3 is devoted to the construction of self-similar functions 
and operators where a self-similar variant of the Laplacian is deduced. This Laplacian gets his physical 
justification in § 4. It is further shown in § 3 that in a continuum approximation this Laplacian takes 
the form of fractional integrals. In § 4 we consider a self-similar quasi-continuous linear chain with 
self-similar harmonic interactions. The equation of motion of this chain takes the form of a self-similar 
wave equation containing the self-similar Laplacian defined in § 3 leading to a dispersion relation having 
the form of the Weierstrass-Mandelbrot function which is a self-similar and for a certain parameter 
range also a fractal function. 

3 Construction of self-similar functions and linear operators 

In this paragraph we define the term "self-similarity" with respect to functions and operators. We call 
a scalar function 0(/i) exact self- similar with respect to variable h if the condition 



is satisfied for all values h > of the scalar variable h. We call (1) the "affine problem"^ where 
is a fixed parameter and A = A^^ represents a continuous set of admissible eigenvalues. The band of 

^ where we restrict here to afRne transformations h' = Nh -\- c with c = 0. 



(f){Nh) = A(j){h) 
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admissible 6 = is to be determined. A function satisfying (1) for a certain N and admissible 
A — represents an unknown "solution" to the affine problem of the form (pN,5ih) and is to be 
determined. 

As we will see below for a given solutions exist only in a certain range of admissible A. 

From the definition of the problem follows that if is a solution of (1) it is also a solution of 
(j){N^h) = A^4>{h) where s e Z is discrete and can take all positive and negative integers including 
zero. We emphasize that non-integer s are not admitted. The discrete set of pairs A*, A^^ are for all 
s G Z related by a power law with the same power S, i.e. A = A^"^ hence A* = (A^*)*^. By replacing A 
and A^ by A^^ and A^"^ in (1) defines the identical problem. Hence we can restrict our considerations 
on fixed values of A^ > 1. 

We can consider the affine problem (1) as the eigenvalue problem for a linear operator A^r with a cer- 
tain given fixed parameter A^ and eigenfunctions (j){h) to be determined which correspond to an admis- 
sible range of eigenvalues A = A^'^ (or equivalently to an admissible range of exponent S — In A/In A"). 
For a function f{x, h) we denote by Af^{h)f{x, h) —: f{x, Nh) when the affine transformation is only 
performed with respect to variable h. 

We assume A, A^ e R for physical reasons without too much loss of generality to be real and 
positive. For our convenience we define the "affine" operator An as follows 

ANf{h) =: f{Nh) (2) 
It is easily verified that the affine operator is linear, i.e. 

An (ci/i(/i) + C2f2{h)) = c.MNh) + C2f2{Nh) (3) 

and 

A%f{h) = f{N'h), s = 0±l,±2,..±oo (4) 

We can define affine operator functions for any smooth function g{T) that can be expanded into a 
Maclaurin series as 

g{T) = f:asT^ (5) 

s=0 

We define an affine operator function in the form 

oo 

g{UN)^E<'se^N (6) 

where ^ denotes a scalar parameter. The operator function which is defined by (6) acts on a function 
f{h) as follows 

oo 

g{aN)f{h) = Y.<^sef{N'h) (7) 

where relation (4) with expansion (6) has been used. The convergence of series (7) has to be verified 
for a function f{h) to be admissible. An explicit representation of the affine operator An can be 
obtained when we write f{h) = f{e^^^) = /(In/i) to arrive at 

AN{h) = e'^'^^TO (8) 
This relation is immediately verified in view of 

AN{h)f{h) = e^°^^/(ei"'^) = /(e'-'^+i-^) = f{Nh) (9) 
With this machinery we are now able to construct self-similar functions and operators. 
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3.1 Construction of self-similar functions 

A self-similar function solving problem (1) is formally given by the series 

oo oo 

m = E ^~'^%fih) = A-^fiN^h) (10) 

s=— oo s=— oo 

for any function f{h) for which the series (10) is uniformly convergent for all h. We introduce the 
self-similar operator 



■N 



E (11) 



that fulfils formally the condition of self-similarity AnT^ = AT/v and hence (10) solves the affine 
problem (1). In view of the symmetry with respect to inversion of the sign of s in (10) and (11) we 
can restrict ourselves to A^ > 1 (A", A e R) without any loss of generality^: Let us look for admissible 
functions f{t) for which (10) is convergent. To this end we have to demand simultaneous convergence 
of the partial sums over positive and negative s. Let us assume that (where we can confine ourselves 
to i > 0) 

lim/(i) = aor (12) 

For t — >• cxD we have to demand that |/(t)| increases not stronger than a power of t, i.e. 

\imf{t)^c^t'' (13) 

t— ►oo 

with ao,Coo denoting constants. Both exponents e IR are allowed to take positive or negative 
values which do not need to be integers. A brief consideration of partial sums yields the following 
requirements for A = A^'^, namely: Summation over s < in (10) requires absolute convergence of a 
geometrical series leading to the condition for its argument AA"~" < 1. That is we have to demand 
S < a. The partial sum over s > requires absolute convergence of a geometrical series leading to the 
condition for its argument A~^A"^ < 1 which corresponds to S > p. Both conditions are simultaneously 
met if 

AT^^ < A = AT'' < A^" (14) 

or equivalently 

Relations (14) and (15) require additionally (3 < a, that is only functions f{t) with the behaviour 
(12) and (13) with [3 < a are admissible in (10). The case ^ — includes for instance certain bounded 
functions \f{t)\ < M such as some periodic functions. 

3.2 A self-similar analogue to the Laplace operator 

In the sprit of (10) and (11) we construct an exactly self-similar function from the second difference 
according to 

(f){x, h) = TNih) {u{x + h)+ u{x - h) - 2u{x)) (16) 



^We also can exclude the trivial case N =1. 
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where «(..) denotes an arbitrary smooth continuous field variable and Tjv(/i) expresses that the affine 

operator A]\[{h) acts only on the dependence on h, that is Aj^{h)v{x,h) = v{x,Nh). We have with 
^ = A"^ the expression 

oo 

(f){x,h)= J2 CMx + N'h)+u{x-N'h)-2u{x)} (17) 

s=—oo 

which is a self-similar function with respect to its dependence on h with A]\f{h)(l){x, h) = Nh) = 
^~^(j){x,h) but a regular function with respect to x. The function (f){x,h) exists if the series (17) is 
convergent. Let us assume that u{x) is a smooth function with a convergent Taylor series for any h. 

Then we have with u{x ± h) — e^'^'^u(x) and u(x + h) + u(x — h) — 2u(x) — (e'^'& + e~^'^ — 2^ u{x) 
which can be written as 

( h d\ (f 

u{x + h) + u{x — h) — 2u{x) = 4sinh^ I o j~ ) ^(•^) ~ TT^^"^) orders /i-^ (18) 

thus a = 2 in criteria (12) is met. If we demand u{x) being Fourier transformable we have as necessary 
condition that 

/oo 
|ti(x)| dx < oo (19) 
-oo 

exists. This is true if \u{t) \ tends to zero as t — > ±oo as \t\^ where (3 < —1. We have then the condition 
that 

/?<0<5 = --i^ <q; = 2 (20) 
InA^ ^ ^ 

We will see below that only (5 > is physically admissible, i.e. compatible with harmonic particle- 
particle interactions which decrease with increasing particle-particle distance. 
The ID Laplacian Ai is defined by 



d^ . . {uix + T) + uix — T) — 2uix)) 

-r^u{x) = hm 

dx^ ' ' r^o r 

Let us now define a self-similar analogue to the ID Laplacian (21) where we put with ^ = 



Aiu{x) = ^u{x) = hm ^ ^— ^ ^ (21) 



A(s,n,t)u(x) —: const lim r (f)(x, r) (22) 

OO 

const liin r"-^ ^ f {u{x + N't) + u{x - N't) - 2u{x)) (23) 



where we have introduced a renormalisation-multiplier with the unknown power A to be determined 
such that the hmiting case is finite. The constant factor const indicates that there is a certain 
arbitrariness in this definition and will be chosen conveniently. Let us consider the limit r ^ by 
the special sequence r„ = N^"'h with n ^ oo and h being constant. Unlike in the ID case (21), the 
result of this limiting process depends crucially on the choice of the sequence r„. We see here that 
the self-similar Laplacian cannot be defined uniquely as in the ID case. We have (by putting in (22) 
const — h^) 

oo 

A^s,N,h)u{x) = Jim N^^'C E («(^ + N'-^'h) + u{x - N'-^'h) - 2u{x)) (24) 

s=—oo 

which assumes by replacing s — n — > s the form 

\5,N,h)u{x) = 0(x, h) lim 7V-(^-^)" (25) 
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which is only finite and nonzero ii X = 6. The "Laplacian" can then be defined simply by 

^is,N,h)u(x) =: jim 7V5"0(x, N'^h) = h) (26) 

or by using (16) and (18) we can simply write^ 

A(,.,.) = Af.ih) sinh^ A j = 4 N-^^ sinh^ A j (27) 

where TN{h) is the self-similar operator defined in (11). The self-similar analogue of Laplace operator 
defined by (27) depends on the parameters S, N, h. We furthermore observe the self-similarity of 
Laplacian (27) with respect to its dependence on h, namely 

^{S,N,Nh) = N^^{S,N,h) (28) 



3.3 Continuum approximation - link to fractional integrals 

For numerical evaluations it may be convenient to utilize a continuum approximation of the self-similar 
Laplacian (27). To this end we put = 1 + e (with < e <^ 1 thus e ^ InN) where e is assumed to be 
"small" and se = v such that dw ~ e and = (1 + e)' ^ . In this approximation iV'^ pa becomes 
a (quasi)-continuous variable when s runs through s e Z. Then we can write (10) in the form 

(j){h) = ^ N-'^f{N'h) ^ - / e-^''f{he'')dv (29) 

s=-oo ^ ■'-'^ 

which can be further written with he^ = t {h > 0) and ^ = dv and t{v — >• —00) = and t{v — >• 

00) = cxD as 



tr^^dr (30) 



-7-1+5 

In this continuum approximation the function 0(/i) obeys the same scaling behaviour as (10), namely 
(j){h\) = \^(j){h) but in contrast to (10) A can assume any continuous positive value. This is due to 
the fact that (30) is holding for iV = 1 -|- e with sufficiently small e > since in this limiting case there 
exists for any continuous value A > an m e Z such that N'^ A. The existence requirement for 
integral (30) leads to the same requirements for f{t) as in (10), namely inequality (15). Application 
of the approximate relation (30) to Laplacian (27) yields 

h' (ii(.r - -) + u{x + -) - 2,/(,t;)) 
^{5,e,h)<^) ~ - | dr (31) 

where this integral exists for /3 < < 5 < 2 and (3 < —1 because the required existence of integral (19) 
and relation (18). By performing two partial integrations and by taking into account the vanishing 
boundary terms at r = and r = 00 for < 5 < 2, we can re-write (31) in the form of a convolution 

of the conventional ID Laplacian ^{x)-, namely 



A^S,e,h)U{x) ^ j^j{\x-T\) — {T)dT (32) 



with the kernel 



^We have to replace 4 — > ^ if the Laplacian acts on a field u{x, t) as in Sec. 4. 



6 



^(ki) = ^^^^kr-', Ml (33) 

where < S < 2 and g{\x\) — — Mn \x\ for 5 = 1. We can further write for 5 7^ 1 (32) in terms of 
fractional integrals 

A^s=2-D,e,H)u{x) ^ ^^^^ ^^^^^ _ (PZ^,, + (-l)^!^"^,) Ai«(x) (34) 

where Aim(x) = -j^u(x) denotes the conventional ID-Laplacian and D = 2 — 5 > which is positive in 
the admissible range of < 5 < 2. For < 5 < 1 the quantity D can be identified with the estimated 
fractal dimension of the fractal dispersion relation of the Laplacian [18] which is deduced in the next 
section. In (34) we have introduced the Riemann-Liouville fractional integral T)'^ which is defined 
by(e.g. [19, 20]) 

K^vix) = ^ Jjx - rf-Mr)dr (35) 

where T{D) denotes the F- function which represents the generalization of the factorial function to 
non- integer D > 0. The F-function is defined as 

/■oo 

r(D) = / T^-^e-^dr , D>0 (36) 
Jo 

For positive integers D > the F-function reproduces the factorial- function r{D) — {D — 1)1 with 
D = 1,2, ..oo. 



4 The physical chain model 

Wc consider an infinitely long quasi-continuous linear chain of identical particles. Any space-point x 
corresponds to a "material point" or particle. The mass density of particles is assumed to be spatially 
homogeneous and equal to one for any space point x. Any particle is associated with one degree 
of freedom which is represented by the displacement field u{x, t) where x is its spatial (Lagrangian) 
coordinate and t indicates time. In this sense we consider a quasi continuous spatial distribution of 
particles. Any particle at space-point x is non-locally connected by harmonic springs of strength ^"^ 
to particles located at x ± N^h, where N > 1 and N e R is not necessarily integer, h > 0, and 
s = 0, ±1, ±2, .. ± oo. The requirement of decreasing spring constants with increasing particle-particle 
distance leads to the requirement that ^ — < 1 [N > 1) i.e. only chains with 6 > are 

physically admissible. In order to get exact self-similarity we avoid the notion of "next-neighbour 
particles" in our chain which would be equivalent to the introduction of an internal length scale (the 
next neighbour distance). To admit particle interactions over arbitrarily close distances N^h 
(s — > — oo, h = const) our chain has to be quasi- continuous. This is the principal difference to the 
discrete chain considered recently by Tarasov [7] which is discrete and not self-similar. 
The Hamiltonian which describes our chain can be written as 

H^- (u'^{x,t) + V{x,t,h))dx (37) 

2 J — OO 



7 



In the spirit of (10) the elastic energy density V{x,t,h) is assumed to be constructed self-similarly, 
namely^ 

V{x, t, h) = ^fW(/i) \{u{x, t) - u{x + h, t)f + {u{x, t) - u{x - h, t)f] (38) 
where TAr(/i) is the self-similar operator (11) with ^ = = to arrive at 

1 oo 

s=— oo 

The elastic energy density V{x, t, h) fulfills the condition of self-similarity with respect to /i, namely 

AN{h)V{x, t, h) = V(x, t, Nh) = r^V(x, h) (40) 

We have to demand in our physical model that the energy is finite, i.e. (39) needs to be convergent 
which yields a = 2 as for the Laplacian (17). To determine (5 we have to demand that u{x,t) is a 
Fourier transformable field^. Thus we have to have an asymptotic behaviour of \u{x ± r, t)| ^ as 
where /5 < — 1 as r — > cxd. From this follows \u{x,t) — u{x ± r, t)p behaves then as \u{x,t)\'^. Hence, 
the elastic energy density (39) is finite if 

0<S<2 (41) 

where P < —1. 

This inequality determines the range of the admissible values of S in order to achieve convergence. 
The equation of motion is obtained from 

^ = (42) 
(where S./Su stands for a functional derivative) to arrive at 



^ - Yl C{Mx,t)-u{x + hN%t)-u{x-hN%t)) (43) 

s=— oo 

= A^s,N,h)u{x,t) (44) 



with the self-similar Laplacian A(j jy./i) of equation (27). As the elastic energy density (39) the equation 
of motion is convergent for 6 being in the interval (41) where /? < —1. We can re-write (44) in the 
compact form of a wave equation 

O{5,N,h)u{x,t)^0 (45) 
where □(5 jv,/i) is the self-similar analogue of the d'Alemhertian wave operator having the form 

^{5,N,h) = \5,N,h) - (46) 

The d'Alembertian (46) with the Laplacian (27) describes the wave propagation in the self-similar 
chain with Hamiltonian (37). The present approach seems to be useful as a point of departure to 
establish a generalized theory of wave propagation in self-similar media. 

Now we determine the dispersion relation, which is constituted by the (negative) eigenvalues of 
the (semi-) negative definite Laplacian (27). To this end we make use of the fact that the displacement 



^The additional factor 1/2 in the elastic energy avoids double counting. 

^This assumption defines the (function) space of eigenmodes and corresponds to infinite body boundary conditions. 
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field u{x,t) is Fourier transformable (guaranteed by choosing (3 < —1) and that the exponentials e^''^ 
are eigenfunctions of the self-similar Laplacian (27). We hence write the Fourier integral 



u{x,t) = ^ r u{k,ty^''dk (47) 

27r J -oQ 



to re-write (44) for the Fourier amplitudes u{k, t) in the form 



and obtain 



{k,t) = -uj\k)u{k,t) (48) 



u;\kh)=A J2 N-'^sm'i-^) (49) 

s=— oo 

The series (49) describes a Weierstrass- Mandelbrot function which is a continuous and for < 
5 < 1 a, nowhere differentiablc function [1, 18]. The Weierstrass-Mandelbrot function (49) fulfills the 
condition of self-similar symmetry, namely 

u^{Nkh) ^ u^{kh) (50) 

where the interval of convergence of the series of the Weierstrass-Mandelbrot function (49) is also given 
by (41). We emphasize that indeed only exponents 5 in the interval (41) are admissible in Hamiltonian 
(37) with the elastic energy density (39) in order to have a "well-posed" problem. 
It was shown by Hardy [18] that for > 1 and ^ = < 1 or equivalently for 

0<6<1 (51) 

the Weierstrass-Mandelbrot function of the form (49) is not only self-similar but also a fractal curve 
of (estimated) non-integer fractal (Hausdorff) dimension D = 2 — S > 1. Figs. 2-4 show dispersion 
curves uo'^{kh) for different decreasing values of admissible < 5 < 1 and increasing fractal dimension 
D. Fig. 1 corresponds to the non-fractal case (5 = 1.2 > 1). The increase of the fractal dimension 
from Figs. 2-4 is indicated by the increasingly irregular harsh behaviour of the curves. In Fig. 4 the 
fractal dimension of the dispersion curve is with D = 1.9 already close to the plane- filling dimension 
2. 

To evaluate (49) approximately it is convenient to replace the series by an integral utilizing a 
similar substitution as in Sec. 3.3 (e ^ InN). By doing so we smoothen the Weierstrass-Mandelbrot 
function (49). It is important to notice that the resulting approximate dispersion relation is hence 
differentiablc and has not any more a fractal dimension D > 1 m the interval (51). For sufficiently 
"small" \k\h {h > 0), i.e. in the long- wave regime we arrive at 

^^^kh) ^ i^Mc (52) 

which is only finite if {\k\hY is in the order of magnitude of e or smaller. This regime which includes 
the long- wave limit A; — > is hence characterized by a power law behaviour cij{k) ~ Const of the 
dispersion relation. The constant C introduced in (52) is given by the integral 

c-2/;(i^d. (53) 

which exists for 5 being within interval (41). 
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This approximation holds for "small" e~lniV7^0(0<e-C 1)^ which corresponds to the limiting 
case that — e" is continuous. In this limiting case we obtain the oscillator density from [11]^ 

Pi.) = (54) 

which is normalized such that p(a;)da; counts the number (per unit length) of normal oscillators having 
frequencies within the interval [00,00 + dco]. We obtain then a power law of the form 



nSh \C 

where S is restricted within interval (41). We observe hence that the power | — 1 is restricted within the 
range 0<| — l<oofor0<5<2, especially with always vanishing oscillator density p{u; — > 0) = 0. 

Wc emphasize that neither is the dependence on k of the Weierstrass-Mandelbrot function (49) 
represented by a continuous '^-dependence nor is this function differentiable with respect to k. 
Application of (54) is hence only justified to be applied to the approximative representation (52) if 
< e <^ 1 thus N = 1 + e is sufficiently close to 1 so that iV* is a quasi-continuous function when 
s runs through s G Z. Hence (54) is not generally applicable to (49) for any arbitrary > 1. We 
can consider (55) as the low-frequency regime a; — > of the oscillator density holding only in the 
quasi-continuous case = 1 -|- e with < e <^ 1. 



5 Conclusions 

We have depicted how self-similar functions and linear operators can be constructed in a simple 
manner by utilizing a certain category of conventional "admissible" functions. This approach enables 
us to construct non-local self-similar analogues to the Laplacian and d'Alembert wave operator. The 
linear self-similar equation of motion describes the propagation of waves in a quasi-continuous linear 
chain with harmonic non-local self-similar particle-interactions. The complexity which comes into play 
by the self-similarity of the non-local interactions is completely captured by the dispersion relations 
which assume the forms of Weierstrass-Mandelbrot functions (49) exhibiting exact self-similarity and 
for certain parameter combinations (relation (51)) fractal features. In a continuum approximation 
the self-similar Laplacian is expressed in terms of fractional integrals (eq. (34)) leading for small k 
(long- wave limit) to a power-law dispersion relation (eq. (52)) and to a power- law oscillator density 
(eq. (55)) in the low-frequency regime. 

The self-similar wave operator (46) with the Laplacian (27) can be generalized to describe wave 
propagation in fractal and self-similar structures which are fractal subspaces embedded in Euclidean 
spaces of 1-3 dimensions. The development of such an approach could be a crucial step towards a better 
understanding of the dynamics in materials with scale hierarchies of internal structures ("multiscale 
materials") which may be idealized as fractal and self-similar materials. 

We hope to inspire further work and collaborations in this direction to develop appropriate ap- 
proaches useful for the modelling of static and dynamic problems in self-similar and fractal structures 
in a wider interdisciplinary context. 

= has to be excluded since it corresponds to iV = 1. 
^The additional prefactor " 2" takes into account the two branches of the dispersion relation (49) (one for fc < and 
one for fc > 0). 
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Figure 2: Dispersion relation uj'^{kh) in arbitrary units for = 1.5 and 5 = 0.7 



